Identification of source terms in wave equation with dynamic boundary conditions

This paper studies an inverse hyperbolic problem for the wave equation with dynamic boundary conditions. It consists of determining some forcing terms from the final overdetermination of the displacement. First, the Fréchet differentiability of the Tikhonov functional is studied, and a gradient formula is obtained via the solution of an associated adjoint problem. Then, the Lipschitz continuity of the gradient is proved. Furthermore, the existence and the uniqueness for the minimization problem are discussed. Finally, some numerical experiments for the reconstruction of an internal wave force are implemented via a conjugate gradient algorithm.     

A higher-order hybrid spline difference method on adaptive mesh for solving singularly perturbed parabolic reaction–diffusion problems with robin-boundary conditions

We propose a hybrid numerical scheme to discretize a class of singularly perturbed parabolic reaction–diffusion problems with robin-boundary conditions on an equidistributed grid. The hybrid difference scheme is developed by using a modified backward difference scheme in time, a combination of the cubic spline and exponential spline difference scheme in space. The proposed scheme uses a cubic spline difference scheme for the discretization of robin-boundary conditions. For the time discretization of the problem, we use the standard uniform mesh while a layer adapted equidistributed grid is generated for the spatial discretization. By equidistributing a curvature-based monitor function, the spatial adaptive grid is able to capture the presence of parabolic boundary layers without using any prior information about the solution. Parameter uniform error estimates are derived to illustrate an optimal convergence of first-order in time and second-order in space for the proposed discretization. The accuracy of the proposed scheme is confirmed by the numerical experiments that underpin the theoretical analysis.     

A non-overlapping domain decomposition dual reciprocity method for solving the forward–backward heat equation in two-dimension

We apply a boundary element dual reciprocity method (DRBEM) to the numerical solution of the forward–backward heat equation in a two-dimensional case. The method is employed for the spatial variable via the fundamental solution of the Laplace equation and the Crank–Nicolson finite difference scheme is utilized to treat the time variable. The physical domain is divided into two non-overlapping subdomains resulting in two standard forward and backward parabolic equations. The subproblems are then treated by the underlying method assuming a virtual boundary in the interface and starting with an initial approximate solution on this boundary followed by updating the solution by an iterative procedure. In addition, we show that the time discrete scheme is unconditionally stable and convergent using the energy method. Furthermore, some computational aspects will be suggested to efficiently deal with the formulation of the proposed method. Finally, two forward–backward problems, for which the exact solution is available, will be numerically solved for two different domains to demonstrate the efficiency of the proposed approach.     

Finite Integral Transform Solutions for Free Vibrations of Rectangular Thin Plates With Mixed Boundary Constraints北大核心CSCD

Analytical solutions, with unique research value, can serve as benchmarks for empirical formulas and numerical methods, a tool for rapid parameter analysis and optimization, and a theoretical basis for experimental designs. Conventional analytical methods, e.g., the Lévy solution method, are only applicable to mechanical problems of plates and shells with opposite simply-supported edges, which, however, may fail to obtain analytical solutions for the issues with complex boundary constraints. In recent years, the finite integral transform method for plate and shell problems was developed to deal with non-Lévy-type plates and shells, but it is still infeasible to solve the mixed boundary constrains-induced complex boundary value problems of higher-order partial differential equations. Herein, for the first time, the finite integral transform method was combined with the sub-domain decomposition technique to solve the free vibrations of rectangular thin plates with mixed boundary constraints. The rectangular plate was first divided into 2 sub-domains according to the mixed boundary constraints, and the 2 sub-domains were solved analytically with the finite integral transform method. Finally, the continuity conditions were introduced to obtain the analytical solution of the original problem. Based on the side spot-welded cantilever plates commonly used in engineering, the free vibration problem of a rectangular thin plate with 1 edge subjected to clamped-simply supported constraints and the other 3 edges free, was analyzed. The obtained natural frequencies and mode shapes are in good agreement with those from the finite element method as well as the solutions in literature, thus verifying the accuracy of the proposed method. The solution procedure of the finite integral transform method can be implemented based on the governing equations without any assumption of the solution form. Therefore, this strict analytical method is widely applicable to complex boundary value problems of higher-order partial differential equations for such mechanical problems of plates and shells. © 2023 Editorial Office of Applied Mathematics and Mechanics. All rights reserved.     

Quantum Bayesian computation

Quantum Bayesian computation is an emerging field that levers the computational gains available from quantum computers. They promise to provide an exponential speed-up in Bayesian computation. Our article adds to the literature in three ways. First, we describe how quantum von Neumann measurement provides quantum versions of popular machine learning algorithms such as Markov chain Monte Carlo and deep learning that are fundamental to Bayesian learning. Second, we describe quantum data encoding methods needed to implement quantum machine learning including the counterparts to traditional feature extraction and kernel embeddings methods. Third, we show how quantum algorithms naturally calculate Bayesian quantities of interest such as posterior distributions and marginal likelihoods. Our goal then is to show how quantum algorithms solve statistical machine learning problems. On the theoretical side, we provide quantum versions of high dimensional regression, Gaussian processes and stochastic gradient descent. On the empirical side, we apply a quantum FFT algorithm to Chicago house price data. Finally, we conclude with directions for future research.     

Absorbing Interface Conditions for the Simulation of Wave Propagation on Non-Uniform Meshes

We proposed absorbing interface conditions for the simulation of linear wave propagation on non-uniform meshes. Based on the superposition principle of second-order linear wave equations, we decompose the interface condition problem into two subproblems around the interface: for the first one the conventional artificial absorbing boundary conditions is applied, while for the second one, the local analytic solutions can be derived. The proposed interface conditions permit a two-way transmission of low-frequency waves across mesh interfaces which can be supported by both coarse and fine meshes, and perform a one-way absorption of high-frequency waves which can only be supported by fine meshes when they travel from fine mesh regions to coarse ones. Numerical examples are presented to illustrate the efficiency of the proposed absorbing interface conditions.     

Asymptotics of Neumann harmonics when a cavity is close to the exterior boundary of the domain

We construct the asymptotics (as ε→0) of solutions to the Neumann problem for the Laplace equation and of the corresponding Dirichlet integral. The problem concerns a three-dimensional domain having two connected components of the boundary at the distance ε>0. To cite this article: G. Cardone et al., C. R. Mecanique 335 (2007).     

具有域内支承的中厚板的边界元法分析

本文借助简化的Reissner理论,利用边界元法对具有域内支承的中厚板进行了分析。建立了相应的边界积分方程,导出了各基本解,讨论了域内支承柱对板的支承情况,编制了计算程序并给出部分算例。     

Interaction of elastic waves with a periodic array of collinear inplane cracks

A boundary integral equation method is applied to the study of the interaction of plane elastic waves with a periodic array of collinear inplane cracks. Numerical results are presented for the dynamic stress intensity factors. The effects of the wave type, wave frequency, wave incidence angle, and crack spacing on the dynamic stress intensity factors are analyzed in detail. The project supported by the Committee of Science and Technology of Shanghai and Tongji University     

基于可靠度的结构优化的序列近似规划算法

基于可靠度的优化的最直观解法是把可靠度和优化的各自算法搭配一起形成嵌套两层次迭代。为改善其收敛性提高计算效率,人们提出了功能测度法、半无限规划法、单层次算法等多种改进方法。本文对传统结构优化界的经典序列近似规划法改造并扩展应用于求解基于可靠度的结构优化问题,构造该问题的序列近似规划模型和求解过程;其核心思想是在每个近似规划子问题中采用近似可靠度指标对设计变量的线性近似,在优化迭代过程中同步更新设计变量和随机空间中的近似验算点坐标,以达到可靠度分析和优化迭代同步收敛的目标。为了算法的实施,还推导出近似可靠度指标的半解析灵敏度计算公式,编制了程序,最终实现与通用软件的连接。论文用算例证实算法的有效性。